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Abstract 

Background: Data on blood flow regulation, renal filtration, and urine output in salt-sensitive Dahl S rats fed on 
high-salt (hypertensive) and low-salt (prehypertensive) diets and salt-resistant Dahl R rats fed on high-salt diets 
were analyzed using a mathematical model of renal blood flow regulation, glomerular filtration, and solute 
transport in a nephron. 

Results: The mechanism of pressure-diuresis and pressure-natriuresis that emerges from simulation of the 
integrated systems is that relatively small increases in glomerular filtration that follow from increases in renal arterial 
pressure cause relatively large increases in urine and sodium output. Furthermore, analysis reveals the minimal 
differences between the experimental cases necessary to explain the observed data. It is determined that 
differences in renal afferent and efferent arterial resistances are able to explain all of the qualitative differences in 
observed flows, filtration rates, and glomerular pressure as well as the differences in the pressure-natriuresis and 
pressure-diuresis relationships in the three groups. The model is able to satisfactorily explain data from all three 
groups without varying parameters associated with glomerular filtration or solute transport in the nephron 
component of the model. 

Conclusions: Thus the differences between the experimental groups are explained solely in terms of difference in 
blood flow regulation. This finding is consistent with the hypothesis that, if a shift in the pressure-natriuresis 
relationship is the primary cause of elevated arterial pressure in the Dahl S rat, then alternation in how renal 
afferent and efferent arterial resistances are regulated represents the primary cause of chronic hypertension in the 
Dahl S rat. 



Background 

Animal models of salt- and/or angiotensin II-induced 
chronic hypertension have revealed shifts in the 
observed pressure-natriuresis and pressure-diuresis rela- 
tionships to higher pressures, as well as altered renal 
blood flow regulation [1-6]. The salt-sensitive Dahl S 
(SS) rat is a widely studied example of an animal that 
develops hypertension, associated with a shift of the 
pressure-natriuresis relationship (relationship between 
sodium excretion and arterial pressure) to higher pres- 
sures, when fed a high-salt diet. When maintained on 
high salt (e.g., 8% NaCl in chow) the kidneys of these 
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animals are found to excrete a given amount of sodium 
per unit time at a higher input arterial pressure than the 
kidneys of control animals fed low-salt diets and of 
strains, such as the salt-resistant Dahl R (SR) rat, that do 
not exhibit salt-induced hypertension. Thus sodium bal- 
ance (dietary sodium input minus sodium excretion) is 
achieved in hypertensive animals at higher pressures 
than in normotensive animals [7]. 

Guyton and Coleman and coworkers hypothesize that 
a shift in the pressure-natriuresis relationship to higher 
pressures is one of the central causal mechanisms of 
chronic hypertension in salt-sensitive hypertension [8]. 
Other investigators suggest that angiotensin II- and salt- 
induced increases in sympathetic nervous activity in the 
vasculature may be a primary causal factor in salt- 
sensitive hypertension while the shift in the renal 
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pressure-natriuresis relationship may not [9-11]. Not 
only is it unclear whether and when the observed 
changes in the pressure-natriuresis relationship are 
causes or consequences of chronic hypertension (or in 
some way both), it remains unclear what specific aspects 
of renal physiology are altered in salt-sensitive hyperten- 
sion, underlying the observed changes in the pressure- 
natriuresis and pressure-diuresis relationships. 

Here we analyze data on blood flow regulation, renal 
filtration, and urine output in SS rats fed on high-salt 
(hypertensive) and low-salt (prehypertensive) diets and 
salt-resistant SR rats fed on high-salt diets. We use a sim- 
ple mathematical model of renal blood flow regulation, 
glomerular filtration, and solute transport in a nephron 
to reveal the minimal differences between the three cases 
necessary to explain the observed data. It is found that 
the differences in renal blood flow, glomerular filtration, 
and pressure-diuresis and pressure-natriuresis relation- 
ships may be explained based solely on differences in 
afferent and efferent arteriole regulation in the hyper- 
tensive (high-salt) SS compared to the salt-resistant SR 
and the low-salt SS controls. 

Sources of data 

Data from the SS and SR rats used for model identifica- 
tion are obtained from Roman [12]. Additional independ- 
ent data from SS and SR rats for model comparison were 
obtained from Roman and Kaldunski [13]. For these data 
sets measurements were made in denervated kidneys 
perfused in vivo with plasma levels of vasopressin, aldos- 
terone, corticosterone, and norepinephrine clamped. 
Data from three experimental groups are analyzed: high- 
salt fed hypertensive SS rats with baseline pressure of 
158 ± 2 mmHg, low-salt fed prehypertensive SS rats with 
baseline pressure of 133 ± 1 mmHg, and high-salt fed SR 
rats with baseline pressure of 124 ± 1 mmHg. Additional 
data for comparison to model predictions are obtained 
from Thompson and Pitts [14] and were obtained in 
normal dogs in which glomerular filtration rate was 
modulated by varying renal arterial pressure. (Data 
from Thompson and Pitts on adrenalectomized and 
sympathectomized dogs show similar trends.) 



B 



25 



co 20 
X 



a> 15 



10 



X 

to 



CD 

rr 



50 



15 



CD 

X 

E 
E 

CD 



CD 

X 



10 



(i-/J>Q, 




- s 

-> 

.,/ 

,-' ✓ S 

y t S 
t 

/ / 

/ J 
s / 


s jX 




SS (high-salt) " 

SS (low-salt) 

SR 



100 150 

P a (mmHg) 



200 



\ 

\ 


















SS (high-salt) 

SS (low-salt) 

SR 






50 100 

Pa 


150 
(mmHg) 


200 



Figure 1 A. Diagram of model. B. Afferent arterial resistances 
governed by Equation (8), with parameter values from Table 1, for 
the three experimental cases analyzed. C. Efferent arterial 
conductivities governed by Equations (6), with parameter values 
from Table 1, for the three experimental cases analyzed. 



Methods 

The mathematical model of renal blood flow, glomerular 
filtration, and mass transport in nephrons (diagrammed in 
Figure 1) is composed of two main components, a model 
for renal blood flow and glomerular filtration and a model 
for mass transport in a representative nephron. Both com- 
ponents are based on modifications made to models pre- 
sented in Chapter 20 of Keener and Sneyd [15]. The blood 
flow and filtration model predicts glomerular filtration 
rate, glomerular pressure, and renal blood flow as func- 
tions of input arterial pressure. The predicted glomerular 



filtration rate and pressure serve as inputs to the nephron 
model, which predicts concentrations of sodium and flows 
in the descending and ascending limbs of the loop of 
Henle and a combined intersitium/ascending vasa recta 
space. Predictions of the overall model are compared to 
data on renal blood flow, glomerular filtration, glomerular 
pressure, efferent capillary pressure, urine flow, and so- 
dium excretion in low-salt fed prehypertensive and high- 
salt fed hypertensive SS and high-salt fed SR rats. 

While the treatment of a single nephron as representa- 
tive of whole kidney is a gross simplification compared to 
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models that capture heterogeneities in loop length and 
the three-dimensional architecture of the tubules and 
vasculature [16-19], the model developed here is appro- 
priate to capture the physiological function analyzed 
here. Models of renal flow regulation and tubuloglomer- 
ular feedback [20-22], transport in the proximal tubule 
and cortex [23-25], medulla [17,18,26-35], collecting duct 
[34,36], and other components [16] have been developed 
to capture much more biophysical detail than the whole- 
kidney model developed here. Yet we are aware that no 
previously developed model, however simplified, com- 
bines renal hemodynamics, filtration, and tubular trans- 
port to simulate and analyze data on whole-kidney 
pressure-natriuresis function. 

Previous models of renal system that capture overall 
kidney function, including the pressure-natriuresis and 
pressure- diuresis phenomena, have been developed 
[37-40]. However, these models do not capture spatially 
distributed transport in kidney, even at the simplified 
level of the model developed here. 

Governing equations for blood flow and filtration 

Flow and filtration along glomerular capillaries is gov- 
erned by the conservation equation for flow in a glom- 
erular capillary, q\(y): 
dqi 
dy 



Kf(P 2 -Pi+n c ),y€(0,L), 



(1) 



where y is the distance along the glomerular capillary, Kf 
is the hydraulic permeability, Yl c (y) is the oncotic pres- 
sure in the plasma, and Pj and P 2 are hydrostatic pres- 
sures in the capillary and Bowman's capsule, respectively. 
This expression assumes that the rate of fluid loss from 
the capillary is linearly proportional to the pressure dif- 
ference driving force and that pressure remains effect- 
ively constant along the length of the glomerular 
capillary. Blood enters the capillary at y = 0 with an input 
oncotic pressure Tl c (y = 0) = IT,- = 28 mmHg and input 
flow *7i(y = 0) = Q 2 . Assuming a linear relationship be- 
tween concentration and oncotic pressure, we have 

n c {y) qi {y) = TliQi. (2) 



Combining Equations (1) and (2) yields 

d ^ = K f (p 2 -p 1+ nfi 

dy \ qi 



(3) 



which is a separable equation that can be solved to yield 
the following relationship between input flow q\ 
(y = 0) = Qi and output flow q i{y = L) = Q e : 



Qi \ 1-a J Qi ■ 



(4) 



where AP = P 1 -P 2 , and a = njAP. The filtrate flow 
(glomerular filtration rate) is computed from the 



difference between input and output blood flows, 

Qf= Qi-Qe- 

Blood flow into the glomerulus satisfies the Ohm's 
Law relationship 



QiRa, 



(5) 



where P a is the input arterial pressure and R a is the affer- 
ent arterial resistance, which is phenomenologically mod- 
eled using the following increasing function of filtration 

where 7?°, R™ m ,n b Q a0 , and C TGF are adjustable para- 
meters, and c a (0) is the sodium concentration in the 
ascending limb at the location where it feeds into the dis- 
tal tubule. The factor Qi"' / {Qi"' + Q io ni ) increases 
smoothly with increasing blood flow, representing an 
autoregulatory vasoconstriction. The factor c a {0)/{C TGF + 
c a {0)) is employed to account for tubular-glomerular feed- 
back: increasing salt concentration in the distal tubule sti- 
mulates vasoconstriction. (The concentration c a (0) is 
obtained from the transport component of the model, 
detailed below.) 

Similarly blood flow out of the glomerulus satisfies the 
Ohm's Law relationship 



Pl-P e =QeRe 



(7) 



where P e is the efferent capillary pressure and R e is the 
efferent arterial resistance, which is phenomenologically 
modeled using the following decreasing function of Pi, 
the input pressure into the arteriole. 



Re(Pl 



i? e max , Pi<(b- R^)/m 
b-mP u P x > {b-R^)/r 



(8) 



where i?™ ax , b, and m are adjustable parameters. Thus 
the efferent arteriole is assumed to contribute to the de- 
creasing behavior of the resistance in direct response to 
increases in pressure beyond a certain cutoff value of Pj. 
Equations (6) and (8) predict that afferent resistance 
increases and efferent resistance decreases as renal per- 
fusion pressure is increased, as illustrated in Figure 1. 
The tubular-glomerular feedback component of the 
model acts in the same direction as the autoregulatory 
factor in Equation (6). As pressure increases, filtration 
rate increases, leading to higher concentrations in the 
distal tubule, decreasing afferent conductivity. 

Roman [12] reported measurements of pressures in 
peritubular capillaries (capillaries downstream of outer 
cortical glomeruli). As the model described here does 
not distinguish between corticomeduUary and juxtame- 
dullary glomeruli, the reported peritubular capillary 
pressures are compared to the model variable P e , effer- 
ent capillary pressure. Data on P e as a function of 
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arterial pressure are used to fit representative function 
for P e (P a ): 



(9) 



which invokes four additional adjustable parameters, P e0 , 
Pel, Pe2> and n pe . 
Filtrate flow satisfies the Ohm's Law relationship 



P 2 -P d = P-dQf, 



(10) 



where Pa is the distal tubule pressure and Rd is the re- 
sistance associated with this pressure drop, assumed 
constant. In the absence of data on distal tubule hydro- 
static pressure, we assume a simple linear proportional- 
ity between arterial pressure and P^: 



Pd = dpdPa, 



(11) 



where a w is set to 0.02, which gives a value of distal tu- 
bule pressure of 2.0-3.6 mmHg over a range of renal 
perfusion pressure of 100 to 180 mmHg. 

Equations (4), (5), (7), and (10), invoking 14 adjustable 
parameters (see Table 1), are solved for the four 
unknowns Q„ Q e , P lt and P 2 to provide model predic- 
tions of these flows and pressures, as well as functions of 
input pressure P a . 

Table 1 Adjustable parameter values 



Governing equations for nephron 

Mass transport in nephrons is represented using a one- 
dimensionally distributed model accounting for flows 
and concentrations in a single representative nephron. 
Thus three-dimensional interactions and the anatomical 
heterogeneity of loop lengths are not taken into account. 
Nevertheless, the model is able to effectively match 
observed pressure-diuresis and pressure-natriuresis rela- 
tionships. The nephron model, diagramed in Figure 1, 
simulates flow and sodium concentration in four regions: 
the descending and ascending limbs of the loop of Henle, 
the collecting duct, and a combined ascending vasa 
recta/interstitium region. Fluid flows in these regions are 
denoted qd, q a , q a and q s ; sodium concentrations are 
denoted C& c a , c c , and c s , where subscripts 'd^ 'aj 'c', and V 
indicate descending limb, ascending limb, collecting 
duct, and interstitial space. After passing through the 
proximal tubule, filtrate enters the descending limb at 
spatial position x = 0; the nephron region is defined over 
the spatial domain x&[0,D], where D = 2 mm is the 
length of the segments of the nephron. 

Fluid transport between the interstitium and the des- 
cending limb is assumed to be linearly proportional to 
the combined mechanical and osmotic pressure driving 
force, Pd + II S - P s + 2RT(cd~ c s ), where Pd is the hydro- 
static pressure in the descending limb, IJ S is the osmotic 



Dahl-R 



Dahl-S 



glomerular hydraulic permeability times length 

resistance associated with distal tubule 

afferent arteriole resistance parameter 

minimum afferent arteriole resistance 

afferent autoregulation parameter 

afferent autoregulation parameter 

maximum efferent arteriole resistance 

efferent arteriole resistance parameter 

efferent arteriole resistance parameter 

TGF concentration parameter 

efferent capillary pressure fitting parameter 

efferent capillary pressure fitting parameter 

efferent capillary pressure fitting parameter 

efferent capillary pressure fitting parameter 

sodium permeability of the descending limb 

hydraulic permeability of the descending limb 

hydraulic permeability of the collecting duct 

maximum sodium reabsorption rate in ascending limb 

apparent Michaelis-Menten constant for sodium reabsorption 



K/L (rnl-rnin _1 -g _1 -rnrnHg -1 ) 
R d (min-g-mmHg-mr 1 ) 
Pp a (ml-min^'-g^-mmHg -1 ) 
f£ lin (ml-min~ 1 .g~'-mmHg~ 1 ) 
Q l0 (ml-min _1 .g H ) 
1/ 

fi™ x (ml-min~ 1 .g~ 1 .mmHg~') 

b (ml-min^-g^'-mmHg" 1 ) 

m (ml-min^-g^-mmHg -2 ) 

Ctgf (mM) 

P e0 (mmHg) 

P e , (mmHg) 

P e2 (mmHg) 

n Pe (unitless) 

H d (ml-min^-g^-mm" 1 ) 

K d (ml-min~ -g~ •mmHg~'-miTr 

K c (ml-min^-g^-mmHg^-mm - 

Pmax (ml-mM-min _1 -mm _1 -g _1 ) 

K m (mM) 



high-salt 

0.0886 

74959 

37.7673 

2.8758 

5.5796 

9.5614 

74185 

37.5995 

0.5269 

25.0 

10.4 

15.1 

136.5 

5.93 

7.70x10 3 

8.3889x10" 

1.8777x10" 

29.172 

50.933 



low-salt 

a 
a 
a 

4.88 
4.28 
a 

13.1 
a 

0.577 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 



high-salt 

a 
a 
a 

9.84 

b 

a 

9.10 

a 

0.54 
a 

8.32 

a 

a 

a 

a 

a 

a 

a 

a 



(a) value same as Dahl-R value. 

(b) value same as Dahl-S value. 



Beard and Mescam BMC Physiology 2012, 12:6 
http://www.biomedcentral.eom/1472-6793/12/6 



Page 5 of 14 



pressure in the interstitium, P s is the hydrostatic pres- 
sure in the interstitium, and c d and c s are the Na + con- 
centrations in the descending limb and the interstitium. 
The factor of 2 multiplying the concentration gradient 
term arises because it is assumed that chloride concen- 
tration equals sodium concentration, and sodium plus 
chloride represent the major contributor to the osmotic 
gradient. With the hydraulic permeability constant K d 
mass conservation yields the equation for q d , the flow in 
the descending limb: 
dq d 



dx 



K d {-AP d + 2RT{c d - c,)),x e [0,D] 



(12) 



where AP d = P d + IJ S - P s , and interstitial osmotic and 
hydrostatic pressures are set to 77 s = 17 mmHg and 
P s = 3 mmHg. 

The ascending limb is assumed impermeable to water, 
and thus flow in the ascending limb, q a , is constant: 

dq a 
dx 

The governing equation for q c , the flow in the collect- 
ing duct is analogous to the equation for q d . 

dq c 



0. 



(13) 



dx 



K c (-AP c + 2RT(c c -c s )) 7 



(14) 



where AP C = P c + TI S - P s and K c is the hydraulic perme- 
ability in the collecting duct. The hydrostatic pressure in 
the collecting duct is assumed to be 1 mmHg lower than 
that in the distal tubule: P c = P d -l mmHg. 
Since total volume is conserved 



(qd 



q a + q c 



(15) 



dq s _ 
dx dx 

Sodium transport is assumed to be governed by passive 
permeation in the descending limb and collecting duct 
and by active transport in the ascending limb. The govern- 
ing equations for Na + flux in descending limb is given by 

d(qdCd) 



dx 



H d (c s - c d ) 



(16) 



where H d is the descending limb permeability. The trans- 
port rate in the ascending limb is given by 



d(q a c a ) 
dx 



~P(Ca) 



(C Y 

y^a, max J 



a +K m ( Ca f + (c a , max y 



(17) 



where the factor 

PmaxCa 



models a saturable process, with P max and K m adjustable 
parameters. The factor 

(C "I 5 

\^a, max J 



(pa) ~t~ (C fl)max ) 



is applied so that the transport rate goes to zero when 
concentrations in the nephron exceed an upper limit. 
Without this factor, concentrations become unbounded 
when the flow in the collecting duct approaches zero. 
Physically, this is because the predicted concentration gra- 
dient increases as flow through the loop of Henle 
decreases. Without this factor the solution becomes math- 
ematically unbounded when pressure drops low enough 
that all of the filtrate is reabsorbed because in this limit q a 
and dqjdx both approach zero, and the only way for d 
{q a c a )/dx to approach a constant value would be for c a 
and/or its gradient to become unbounded. Since the con- 
centration gradient drives fluid loss from the descending 
limb and the collecting duct, increases in the concentra- 
tion gradient lead to further decreases in flow. The Hill 
coefficient of 5 in this multiplying factor is also arbitrarily 
assigned so that transport rapidly approaches zero when 
c a exceeds C a max . The value of the fixed parameter C aimax 
set to 500 mM, so that the maximal Na + concentration 
achieved at low flows is approximately 800 mM, associated 
with an approximately 5-fold magnification of the input 
concentration of c d (0) = 150 mM. (See below.) For pres- 
sures and flows that result in urine flows that are approxi- 
mately equal to and greater than the baseline values, c a 
remains well below C ajmax and the behavior of the model 
is not sensitive to the values of these fixed parameters. 

Sodium reabsorption in the collecting duct is not ex- 
plicitly accounted for in the model and the equation for 
Na + flux in the collecting duct is 



dx 



0. 



(18) 



This simplifying assumption is discussed below. Salt 
transport in the interstitial space combines active trans- 
port and passive permeation processes: 



dx 



+P(c a ) - H d (c s - c d ) - H c (c s - c c 



(19) 



Equation (19) assumes that the combined interstitial and 
vasa recta space gathers the sum of the fluxes from the 
other structures. Thus, as expressed by Keener and Sneyd, 



dx 



dx 



(qdCd + q a c a + q c c c 



(20) 



The boundary conditions for input into the descending 
limb assume that input concentration is equivalent to 
plasma sodium concentration of 150 mM and input flow 
is proportional to Qp the glomerular filtration rate: 



q d (0)=PQ f , 
c d (0) = 150 mM, 



(21) 



where Qf= Q, - Q e is determined as a function of arterial 
pressure by the renal blood flow and filtration model 
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component and yS = 0.33 is a fixed constant accounting 
for reabsorption by the proximal tubule. Thus, the 
model assumes constant glomerulotubular balance and 
isotonic reabsorption of water and sodium from the 
proximal tubule [41]. 

The boundary conditions for the ascending limb are 
obtained from the assumption of continuity of concen- 
tration and flow at the turn of the loop of Henle: 



q a (D) = -q d (D), 
c a (D) = c d (D). 



(22) 



Similarly, the ascending limb feeds into the collecting 
duct 



1c(0) 
Q(0) : 



,(o). 



The interstitial flow boundary condition is [15] 

q s {D) = 0 

q d (0)+q s (0) =q c (D). 



(23) 



(24) 



Equations (20), (21), (22), and (23) provide eight 
boundary conditions for the eight first-order differential 
equations described above. The eighth condition comes 
from conservation of total fluid flow [15], requiring that 
flow into the system at x = 0 equal flow out of the system 
at x = D. 

Numerical discretization of the nephron model is 
described in the Appendix. 

Since sodium reabsorption occurs in the model only in 
the ascending limb and the outflow of the ascending 
limb feeds directly into the collecting duct, the model 
does not explicitly account for sodium reabsorption in 
the distal tubule or the collecting duct. Thus all sodium 
reabsorption processes are represented by the ascending 
limb sodium transport rate P(C a ). This simplifying ap- 
proximation is justified by the fact that during formation 
of either concentrated or dilute urine, the majority of so- 
dium reabsorption occurs via the ascending limb. Lump- 
ing all reabsorption processes into Equation (17) helps 
keep the model tractable and identifiable. Adding add- 
itional processes would add additional uncertainty in 
parameter values that would not be justified given the 
available data or yield any additional insight into the op- 
eration of the integrated model. 

Results 

Model identification 

Predictions of the renal flow and filtration model com- 
ponent are compared in Figure 2 to data on blood flow, 
filtration rate, glomerular pressure, and efferent pres- 
sure, as functions of arterial pressure in the SS (high-salt 
and low-salt) and SR (high-salt) rats. Predictions of urine 
output (Q u = Q c {x = D)) and sodium excretion (Q c {x = D). 



C c (x = D)) are compared in Figure 3 to data from SS 
(high-salt and low-salt) and SR (high-salt) rats as func- 
tions of arterial pressure. Data plotted in both Figures 2 
and 3 used for model identification are obtained from 
Roman [12]. Additional independent data on urine out- 
put and sodium excretion in the SR rat were obtained 
from Roman and Kaldunski [13]. 

The 19 adjustable parameters invoked in this model 
are not identifiable for a given experimental group based 
on the six data sets (renal blood flow, filtration, glomeru- 
lar pressure, efferent pressure, urine output, and sodium 
excretion as function of renal perfusion pressure) repre- 
sented in these figures. However, the combined data set 
of pressures and flows versus arterial pressure for three 
experimental cases — prehypertensive (low-salt fed) and 
hypertensive (high-salt) SS rats and salt-resistant (high- 
salt) SR rats — provides independent data that can be 
compared to sixteen model-predicted functions of P a : Q„ 
Qfi Pi, P e > Qui and Q u . C c {x = D) under both high- and 
low-salt conditions for the SS and high-salt conditions 
for the SR. (Data onP ; and P e as functions of P a for the 
SS on low salt are not available; data for the SR were used 
to parameterize Equation (9) to represent P e (P a ) for the 
SS on low salt.) If we assume that most model para- 
meters attain the same values for all three groups, it is 
possible to determine identifiable parameter sets and to 
determine the minimal set of differences between the 
two conditions that is able to explain the observed data. 

Specifically, if it is assumed that only two afferent arter- 
ial flow regulation parameters Q° and R° a and the efferent 
arterial flow regulation parameters m and 7?™ ax are differ- 
ent between the SR and SS (low-salt) cases and that only 
the parameters R° a , RT^i vn, and P e0 are different between 
the SS low-salt and high-salt cases, then there are a total 
of 27 adjustable parameters that can be estimated by 
matching data to the 16 model-predicted functions in 
Figures 2 and 3. These parameters that are allowed to at- 
tain different values between the experimental groups 
govern how afferent arterial resistances are regulated in 
the model and do not directly affect glomerular filtration 
or transport in the nephron. Model simulations asso- 
ciated with the parameter values listed in Table 1 are 
plotted in Figures 2 and 3. The data on flows and pres- 
sures in SS and SR rats shown in Figure 2 are effectively 
captured by the model, with the exception of glomerular 
filtration rate in the SR rat. (The apparent mismatch be- 
tween model prediction and reported data on glomerular 
filtration rate in the SR rat is discussed below.) 

The predicted trends in afferent conductivity shown in 
Figure IB may be compared to the measurements of 
Takenaka et al. [42], who observed that: (1.) afferent 
arterioles from low-salt SS animals maintain higher dia- 
meters at low pressures than those from high-salt ani- 
mals; and (2.) arterioles from low-salt SS animals show a 
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Figure 2 Model predictions of renal blood flow and filtration compared to data from Dahl S (SS) rats on high-salt and low-salt diets 
and Dahl R (SR) rats on high-salt diets. Data on renal flow, filtration rate, glomerular pressure, and pertibular efferent pressure are plotted as 
circles for SS high-salt fed animals, triangles for SS low-salt fed animals, and filled squares for SR animals. Model predictions based on the 
parameter values from Table 1 are compared to observed data. The difference between the high-salt and low-salt case is captured primarily by 
differences in afferent arterial resistance control parameters. Data in panel A are obtained from Figure 1 of Roman [12]; panel B from Figure five 
of Roman [12]; panel C and D from Figure six of Roman [12]. 



stronger constriction in response to increasing pressure 
than those from high-salt animals. These results are con- 
sistent with our model predictions. However, the obser- 
vations of Takenaka et al., which made use of an isolated 
buffer perfused hydrophonetic kidney preparation, 



demonstrated total abolishment of the autoregulatory re- 
sponse in afferent arterioles in the high-salt case. An 
exact match between the model and the in vitro data of 
Takenaka et al. is not expected because the nature of the 
experiments of Takenaka et al. altered any sheer- 
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Figure 3 Predicted pressure-natriuresis and pressure-diuresis relationships. Urine output (Q u = Q c (x = D)) and sodium excretion (Q c (x = D). 
C c (x=D)) are plotted as functions of arterial pressure, for the Dahl salt-sensitive group on high-salt diet (circles) and low-salt diet (triangles) and 
for the Dahl salt-resistant group as filled symbols (squares and diamonds). Data for urine output and sodium excrection are obtained from 
Figure 3 of Roman [12] and Figure 5 of Roman and Kaldunski [13]. The data for the Dahl R group plotted as diamonds are obtained from Roman 
and Kaldunski [13]; all other data are obtained from Roman [12]. Model predictions for all cases use parameter values defined in Table 1. 
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dependent component of physiological diameter regula- 
tion and abolished tubular-glomerular feedback, and be- 
cause the data of Roman show a clear, if blunted, 
autoregulation of renal blood flow and filtration in the 
high-salt animals. 

The acute changes in sodium excretion and urine out- 
put in response to changes in renal perfusion pressure 
plotted in Figure 3 are termed the pressure-natriuresis 
and pressure-diuresis curves. These acute responses 
should not be confused with long-term relationships be- 
tween pressure and sodium excretion and urine output, 
which are influenced by a number of hormonal, neural, 
and remodeling processes not accounted for here. Here, 
the acute pressure-natriuresis and diuresis phenomena 
are effectively reproduced by the model. Since the glom- 
erular filtration and nephron transport parameters are 
held fixed for all experimental groups, the differences in 
afferent and efferent arteriole tone are responsible for 
greatly diminished rates of urine output and increased 
rate of sodium reabsorption in the SS (on high and low 
salt) compared to the SR. 

Predictions of concentration and flow profiles in the 
nephron, based on the nephron model, are illustrated in 



Figure 4. The upper panel plots model predictions asso- 
ciated with an arterial pressure of 125 mmHg, near the 
baseline pressure of 126 ± 1 mmHg observed in the 
(high-salt) SR rats [12]. The lower panel plots model 
predictions associated with lowering the input pressure 
to 95 mmHg. Although the differences in input pressure 
and flow between the upper and lower panels are small, 
the predicted model behaviors show a major qualitative 
difference. The slightly lower input flow for the lower 
pressure simulation results in collecting duct flow that 
drops to near zero at the outlet at x = 2 mm. Also, at 
lower flow the concentration gradient is greater than at 
the higher flow. At arterial pressures 95 mmHg and 
below, the maximal concentrations at x = 2 mm are ap- 
proximately 500 mM, over a three-fold increase of the 
input concentration of 150 mM. 

To summarize the findings of comparing model pre- 
dictions to data from Roman [12] on high-salt SR and 
hypertensive (high-salt) and prehypertensive (low-salt) 
SS rats, the observed differences in renal function may 
be explained primarily by differences in the control of 
afferent and efferent resistance and in sodium reabsorp- 
tion kinetics. 
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Figure 4 Model-predicted sodium concentration and flow profiles in the nephron model. Simulations are conducted using the parameter 
values for the Dahl R rat (Table 1). A. & C. Sodium concentrations as functions of distance along the nephron are plotted for the descending and 
ascending limb, collecting duct, and interstitium. B. & D. Flows as functions of distance along the nephron are plotted for the descending and 
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pressure the concentration gradient steepens and output flow drops to near zero. 



Beard and Mescam BMC Physiology 2012, 12:6 
http://www.biomedcentral.eom/1472-6793/12/6 



Page 9 of 14 



Discussion 

Mechanisms of pressure-natriuresis and pressure-diuresis 

Using a simple mathematical model to simulate blood 
flow regulation, glomerular filtration, and medullary sol- 
ute transport in the kidney, we have analyzed data from 
Dahl S and Dahl R rats to investigate the potential 
mechanistic underpinnings of renal function observed in 
these animals. While the pressure-natriuresis and 
pressure-diuresis relationships illustrated in Figure 3 
have long been recognized as playing a central role in 
the long-term control of blood pressure [43], the bio- 
physical mechanisms underlying these phenomena have 
not been fully resolved. One school of thought maintains 
that because renal blood flow and glomerular filtration 
rate do not change over a wide range of arterial pressure, 
the observed decrease in sodium and water reabsorption 
associated with an increase in pressure could not be sub- 
stantially impacted by an increase in the rate of filtrate 
delivery to nephrons. Observations of reduced sodium- 
hydrogen exchanger activity in the proximal tubule [44] 
and increased medullary blood flow [45,46] in response 
to acute and chronic increases in blood flow point to 
mechanisms that reduce sodium reabsorption. Another 
widely held view is that both acute and chronic increases 
in arterial pressure increase filtration rate and thus dir- 
ectly increase sodium and water excretion through sim- 
ple mechanical transduction. In a textbook explanation 
of the pressure-diuresis phenomenon, this primary 
mechanical response is enhanced by other contributing 
mechanisms, including the renin-angiotensin system, 
changes in medullary blood flow, regulation of proximal 
tubule sodium transport [41]. 

Our model simulations, as well as the data analyzed 
here, are consistent with the pressure-natriuresis and 
pressure-diuresis phenomena emerging from the mech- 
anical relationships between renal pressure, flow, and fil- 
tration. Specifically, in the model increasing arterial 
pressure causes increased glomerular pressure, which 
causes increasing filtration rate. For the SS rat data sets, 
an increase in glomerular pressure of 20-30% over the 
observed pressure range results in an increase in filtra- 
tion rate of 30% in the low-salt case and almost 90% in 
the high-salt case. When pressure increases from 100 
to 180 mmHg, filtrate flow increases from 700 to 
over 900 ul.min _1 .g _1 while urine output increases 
from 10 to 60 ul.min~ 1 .g~ 1 in low-salt case. In the 
high salt case, filtrate flow increases from approximately 
480 to 900 ul.min -1 while urine output increases from 
6 to 68 ul.min ^g 1 over the pressure range of 120 to 
200 mmHg. Thus the slope of filtrate flow (Qf) versus 
arterial pressure can be substantially steeper than the 
slope of urine output (Q u ) versus arterial pressure, even 
over the pressure range for which blood flow is autore- 
gulated. For these cases the relative change in urine 



output over the pressure range is much greater than 
the relative change in filtrate flow because at the lowest 
pressures nearly all of the filtrate is reabsorbed. 

In contrast, the SR data show relatively little increase 
in filtration over the observed pressure range for the 
three data points in Figure 2B. For this case, an increase 
in filtration of approximately 50 ul.min ^g 1 is asso- 
ciated with an increase of 65-80 ul.min _1 .g _1 in urine 
output. For this case the model is not able to capture 
the nearly constant Qf as a function of P a because the 
glomerular pressure is observed to increase from 44 to 
53 mmHg over the same arterial pressure range. Recall 
that the driving force for filtration is hydrostatic pressure 
difference minus the oncotic pressure of approximately 
28 mmHg. Since the 8 mmHg increase in glomerular 
pressure over the observed range of arterial pressure 
represents an approximately 30% increase in driving 
force for filtration, the model tends to under-represent 
the slope of Pj versus P a while over representing the 
slope of Qf versus P a . It is unclear how to resolve the 
substantial differences in driving force for filtration with 
the apparently constant filtration rate observed in the SR 
rat. The model predicts that Qf increases roughly 20% 
over the observed 50 mmHg range of arterial pressure, 
while measurements in the SS rat and in other rat 
strains and other species show increases of anywhere 
from 10% to greater than 20% over the pressure range of 
autoregulated blood flow [47-49]. 

The relationship between sodium excretion and glom- 
erular filtrate rate is further explored in Figure 5 by 
comparing the model predictions of these variables to 
the data of Thompson and Pitts [14]. Here, model pre- 
dictions and data are plotted as percent of control since 
the data are obtained from dog and the model is para- 
meterized for the SR rat. Note that this comparison 
represents a model prediction where no parameter ad- 
justment has been done to match the data. The non- 
linear nature of the relationship is effectively captured 
by the model, where relatively small increases in filtra- 
tion rate can effect relatively large changes in sodium ex- 
cretion. Furthermore, the simplified model reveals the 
extent to which mechanisms not included in the model 
may be important contributors to the physiological phe- 
nomena explored. Specifically, Figure 5 shows that a 
20% increase in filtration rate above baseline level elicits 
a 100% increase in model-predicted urine output rate, 
for the SR parameter set. 

The mechanistic explanation for the pressure-diuresis 
and pressure-natriuresis phenomena that emerges is illu- 
strated in Figure 6. The upper panel plots conceptualized 
curves representing glomerular filtration flow and urine 
flow as functions of arterial pressure. Consistent with the 
available data, the slopes of glomerular filtration flow and 
urine flow versus pressure are of the same order of 
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Figure 5 Relationship between urine output and sodium 
excretion and glomerular filtration rate. A Relative change in 
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the arterial pressure range studied in Figure 2. B. Relative change in 
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Data from Thompson and Pitts [14] from normal dogs are plotted 
for comparison. All calculations are for the Dahl R rat parameter set. 



magnitude over the autoregulated range, here taken to be 
P a = 100 to 160 mmHg. If the slopes are the same over 
this pressure range, then Q u can be approximated as Qf 
minus a constant reabsorbed volume. Assuming that so- 
dium concentration remains approximately constant at 
arterial pressure above the baseline 100 mmHg, the 
pressure-diuresis relationship of the bottom panel is 
obtained. While this conceptual model is highly simpli- 
fied, it does effectively illustrate the basic mechanism 
that emerges from our mathematical model: since glom- 
erular filtration flow is much larger than urine flow, a 
relatively small increase in glomerular filtration can cause 
a relatively large increase in urine output. Thus, this ex- 
planation requires that glomerular filtration does in- 
crease, albeit slightly, as renal arterial pressure is acutely 
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natriuresis. Idealized curves are used to illustrate the hypothesized 
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(QJI, and sodium excretion following acute changes in renal arterial 
pressure (P a ). 



increased. If, as has been hypothesized, glomerular filtra- 
tion remains exactly constant as arterial pressure is 
acutely increased, then this mechanism cannot explain 
the observed pressure-diuresis and pressure-natriuresis 
relationships. 

The data from Thompson and Pitts, as well as the data 
of Roman analyzed in Figures 2 and 3, indicate that the 
slope of Qf versus P a may be lower in normal animals 
than that captured by the model and that additional sec- 
ondary mechanisms may be necessary to satisfactorily 
explain the pressure-diuresis/natriuresis phenomenon. 
Clearly, mechanical transduction is not the only mech- 
anism at work in vivo. Yet even without a model simula- 
tion, it is apparent from the raw data that a given 
change in pressure can induce a greater change in glom- 
erular filtration than in urine output over the pressure 
range for which blood flow is autoregulated. The model 
reveals the extent to which the relationship between 
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acute changes in arterial pressure and glomerular filtra- 
tion can explain the observed pressure-natriuresis and 
pressure-diuresis relationships. These findings highlight 
the direct effects of pressure on influencing urine pro- 
duction by delivering increased filtrate to the proximal 
tubule. The ability of the model to match observed rela- 
tionships between arterial pressure and glomerular pres- 
sure, urine flow, and sodium excretion depends on the 
predicted increase in glomerular filtration with pressure 
that is not apparent in the SR data set. For the model to 
capture the phenomenon of constant glomerular filtra- 
tion over the arterial pressure range of 100 to 150 mmHg 
observed in the SR group would require the introduction 
of some (unknown) mechanism that reduces glomerular 
hydraulic permeability in response to acute increases in 
pressure. Furthermore, since it does not account for hor- 
monal or nervous factors, or changes in medullary blood 
flow and transporter activities, the model reveals that 
these factors are not necessary to explain the acute 
pressure-natriuresis and pressure-diuresis phenomena, 
at least in the SS rat on low- and high-salt diets. 

Physiological differences between SS and SR groups 

In addition to revealing insight into how sodium excre- 
tion and urine output are influenced by perfusion pres- 
sure, model analysis reveals potential mechanistic 
underpinnings of differences in renal function observed 
in SR and SS rats when fed on low-salt versus high-salt 
diets. Our strategy identifies the minimal differences be- 
tween model parameterizations necessary to explain the 
data from these groups. Specifically, it was found that 
differences in five parameters associated with blood flow 
control (see Table 1) are able to explain a host of differ- 
ences in renal function observed among the three groups 
(see Figures 2 and 3). 

In developing the model presented here and deter- 
mining the difference in parameter values necessary to 
explain the groups, the goal is not to capture all rele- 
vant physiological processes impacting renal function 
and blood pressure regulation in the rat. Indeed sev- 
eral mechanisms important to the renal response to 
changes in blood pressure, including changes in prox- 
imal tubule sodium transport [44] and inner medullary 
blood flow [45], are not accounted for. In contrast, by 
focusing on a well-defined set of identifiable physio- 
logical processes, we are able to determine a minimal 
set of processes to explain the data and what differ- 
ences in those processes are necessary to explain the 
different experimental cases. 

The different parameterizations used to explain the 
different experimental groups point to increases in affer- 
ent resistance and decreases in sodium transport rate as 
one moves from lower-pressure to higher-pressure ani- 
mals. The increase in afferent resistance is able to 



explain all of the qualitative differences between 
observed data on renal function — lower flows, filtration 
rates, and glomerular pressure as well as the shift in the 
pressure-natriuresis and pressure-diuresis relationships 
in higher pressure animals. (Since the data analyzed here 
are obtained from denervated kidneys, the predicted dif- 
ferences in afferent arterial tone cannot be explained 
based on differences in sympathetic tone, unless chronic 
differences in sympathetic tone had the effect of chron- 
ically altering afferent arterial tone in a way that is 
reflected in denervated kidneys.) 

This observed shift (compared to lower pressure 
controls) of the pressure-natriuresis relationship to 
higher pressure necessarily occurs in hypertension. 
This is because net sodium balance, by definition, 
must occur at a higher pressure in hypertension than 
in normotension. The view that the chronic pressure- 
natriuresis relationship (also called the renal function 
curve) observed in normal animals is effectively infin- 
itely steep implies that the kidney can maintain blood 
volume and sodium at nearly constant levels in re- 
sponse to small changes in pressure associated with 
salt-loading and volume expansion [50]. Furthermore, 
the steepness of the renal function curve forms the 
basis of the theory that chronic hypertension in 
angiotensin- and salt-induced models is caused by 
renal dysfunction leading to decreased sodium excre- 
tion at a given arterial pressure [1-6]. 

While it is debated whether and/or when renal dys- 
function represents the primary cause of chronic hyper- 
tension in the SS rat (and other animal models) [51,52], 
it is clear that the high-salt diet does cause a shift in the 
acutely measured pressure-natriuresis relationship in the 
SS rat, as illustrated in Figure 3. These changes are 
shown here to be underpinned by changes in renal affer- 
ent and efferent arterial resistance. If indeed a shift in 
the pressure-natriuresis relationship is the primary cause 
of elevated arterial pressure in the Dahl S rat, then alter- 
nation in how renal afferent and efferent arterial resis- 
tances are regulated represents the primary cause of 
chronic hypertension in the Dahl S rat. 

Assumptions and simplifications of the model 

As discussed in the methods section, the developed 
model for blood flow, glomerular filtration, and so- 
dium transport in the proximal tubule, nephron, and 
collecting duct developed here is relatively unsophisti- 
cated compared to a number of previously developed 
models of the three-dimensional architecture of the 
tubules and vasculature [16-19], models of renal flow 
regulation and tubuloglomerular feedback [20-22], 
transport in the proximal tubule and cortex [23-25], 
medulla [17,18,26-35], collecting duct [34,36], and 
other components [16]. Nor does the model account 
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for different roles for the inner and outer medulla, in 
terms of tubular function or blood flow, or explicitly 
for the transport of urea or other solutes. Yet despite 
the simplifications, the model represents the only avail- 
able computational model of whole-kidney function 
that integrates blood flow regulation, glomerular filtra- 
tion, distributed solute and volume fluxes along a 
nephron, and tubuloglomerular feedback, to compare 
whole-kidney pressure-natriuresis and pressure-diuresis 
phenomena to experimentally observed data. 

To a certain degree the level of simplification adopted 
by the model is justified by the nature of the data ana- 
lyzed here and the specific questions addressed by the 
model analysis. The appropriate level of complexity 
represented by a model is the lowest (most simple) that 
can capture the biophysical processes underlying the 
phenomena studied. Based on this standard, the present 
model may be judged as a reasonable, if imperfect, sim- 
plification. The most obvious feature that the model 
does not capture well is the phenomenon of nearly con- 
stant filtration rate observed in the SR group. It is not 
known what anatomical and physiological features not 
represented in the current model are critical to improve 
the behavior of the model in comparison to this observa- 
tion. While the model makes a number of simplifying 
assumptions, it is not clear that relaxing any one of 
those simplifying assumptions would explain the appar- 
ent disconnect between driving force for filtration and 
filtration rate in the SR group. What is clear is that the 
basic model introduced here represents a useful frame- 
work for exploring such questions in a systematic 
matter. 



Conclusions 

Analysis of data on renal blood flow, filtration, pressure- 
diuresis and pressure-natriuresis phenomena in Dahl S 
and Dahl R rats using a simple mathematical model 
reveals a hypothetical mechanistic explanation for the 
observed pressure-diuresis and pressure-natriuresis rela- 
tionships. Idealized curves plotted in Figure 6 illustrate 
the hypothesized relationships between glomerular fil- 
tration flow (Qf), urine output (Q M ), and sodium excre- 
tion following acute changes in renal arterial pressure 
(P a ). Increasing pressure is associated with a relatively 
small increase in glomerular filtration, which increases 
delivery of filtrate to the nephron, leading to increased 
urine production. This simplified conceptual model 
requires that glomerular filtration increases slightly as 
renal arterial pressure is acutely increased. Furthermore, 
differences between Dahl salt-sensitive (SS) and salt- 
resistant (SR) rats in renal filtration and urine produc- 
tion are explained in terms of difference in blood flow 
regulation. 



Appendix: Discretization of nephron equations 

Equations (12), (13), and (14) are discretized using finite 
differences 



7=1 



Qa = ~Qd N> 



Q <i> = q <n> + AxKc (-AP e +2RT(C^ > -C^>)) 
/•=i 

Q <,> = Q <*> _ ( Q<«> + Q <i> + Q<i>^ (25) 

where z'= 1, 2, . . .,N is the element index and Cfi = <7d(0) 
is the input flow into the descending limb. The discrete 
variables Qd l> , Qa^i Qc l> > and Q^ l> are numerical approx- 
imations for the continuous variables qd(x), q a {x), q c {x), 
and q s (x). The equation for Q^ l> is based on mass con- 
servation and the boundary condition > = 0. 

The concentrations satisfy numerical approximations 
of Equations (16) for the descending limb: 

QTCr - Q° d C° d = AxH d {Cf> - Cf>),i = 1 

Q <(>/-"■<!> /-i<j-l>/^<i-l> 
d W Vrf w 

= AxH d (Cf> -C< i> ),i=2,3,...,N (26) 
Equation (17) for the ascending limb: 

C<'> = C<">+**P(C< t >),i = N 
Ax 



C< 1> = C< i+1 > + -^P{C<»), i = 1, 2, . . . , N - 1 



(27) 



Equation (18) or the collecting duct: 



Qf>Cf> - {-Q a )C< x > = AxH c (Cf> - C<<>) ,i=l 

= AxH c (<:<'> - Cf l> ) , i = 2, 3, . . . , N (28) 

and Equation (18) for the interstitium: 

Q0 = Q<1>-[Q0_(Q<1> + Q fl + Q< 1 >)] 
Q? >c ? +1> ~ Q° s Cr = AxH c {Cf> - c<<>) 

+ AxH d (qf> - C s <!> ) + AxPiC^) ,i=l 
Qf> C <i+i> _ Qf-^C^ = AxH c (C< l> - C< 1> ) 

+ AxH d {Cf > -Cf > )+AxP{Cf > ), i=2, 3,. . .,N-1 
Qf>Cf - Qr'>C<'> = AxH c (Cf> - C s «>) 

+ AxH d (qp - Cf>) + AxPiC?*) ,i = N (29) 

The input concentration and flow are obtained from 
the boundary conditions Q d = qd(0) and Cj = Cd(0). Qs <N> 
and Cf are the input vasa recta flow and concentration 
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(at x = D); Qs = q s (D) = 0. As long as Qs = 0, the value 
of Cf is arbitrary. 

These equations are solved using an iterative method. 
Computer codes for the model can be obtained by con- 
tacting the author. 
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